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This paper completes the work, initiated in 00, further 
referred as Parts I and II, concerning to Dirac's quantiza- 
tion of Nambu-Goto theory of open string, formulated in 
the space-time of dimension d = 4. Here we perform more 
detailed study of Gribov's copies in the classical mechan- 
ics and determine the quantum spectrum of masses for the 
arbitrary spin case. 

1 Gribov's copies 

Gribov's copies are multiple intersections of the orbit of 
gauge group and the surface of gauge fixing condition. In 
Part II we have shown that the considered version of string 
theory possesses Gribov's copies, related with the singular 
points of a specially constructed vector field on the sphere. 
Therefore, the following topological invariants can be intro- 
duced for Gribov's copies. 

Definition 1: let the phase space M be a smooth orientable 
manifold. Let the orbit of gauge group G and the surface 
of gauge fixing condition F be its smooth orientable sub- 
manifolds with dim(G) = codim(F). Let P be the point 
of their transversal intersection. This means that in point 
P the tangent spaces to G and F span the tangent space 
to M. Let f(F,P), f(G,P) and f(M,P) be the bases in 
the tangent spaces, defining the orientation of F, G and M, 
evaluated in point P. The index of intersection of F and G 
in point P is defined as a number v, equal to (+1) if the ba- 
sis (t(F, P) , t(G , P)) has the same orientation as the basis 
f(M,P), and equal to (—1) if the orientations are opposite. 

Definition 2: non-degenerate singular points of 2- 
dimensional vector field are focus, center, node and saddle. 
Standard topological classification introduces the index of 
singular point u, equal to (+1) for the cases (focus, center, 
node) and equal to (—1) for the saddle. 

Lemma 1 (topological charges of Gribov's copies): 
these two definitions coincide. 

The explicit expressions for the vector fields are given by 
Eq.(12) in Part II. There is also an alternative definition: 

■ exp [22* (L(a) - (Q(a) - X)e s )] , (1) 

where X = <j> dL(a)Q(a) / L t defines an average position of 

the curve Q(o~) and coincides with the definition of mean 
coordinate in CMF given in Part I. The difference of Q and 
PII-Eq.(12) consists in ^-dependent phase factor e mv{ - S3) , 

ip{es) = 2n(X — Q(Q))ez/L t . According to the Part II, the 
evolution of a n is the phase rotation. Thus, the phase fac- 
tor actually introduces a difference of "local time" for the 
evolution of vector fields on the sphere. 



Lemma 2 (phase factor): the factor e'" v< - e3 ^ preserves the 
orbits of Gribov's copies and changes their evolution param- 
eter from leg's PII-Eq.(13) to the natural one (length of the 
curve) . 

Further we will use the definition Q , more convenient to 
study the evolution of Gribov's copies. 
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Fig.l. Gribov's copies. These images are created by a com- 
puter program for the visualization of vector fields, courtesy of 
Wasil Urazmetov, IHEP, Protvino, Russia. A paper on the visu- 
alization techniques used in this program will be published else- 
where. 

To investigate the properties of Gribov's copies in more 
detail, we perform their visualization. The vector fields 
of figQ are displayed on a spherical map (longitude, lati- 
tude), where the top and bottom lines correspond to auxil- 
iary poles N',S'. These poles are specially selected in regular 



points of the vector field, to prevent the image distortion in 
the vicinity of the original poles N,S, where the singulari- 
ties are located. The equator E for the poles N,S is shown 
on the images by dashed line. The following singularities 
of the vector field are visible in these images: nodes - as 
spots of black/white color, dependent on the field direction 
(source/sink); saddles - as cross-shaped black- white spots 
with a hyperbolic structure of the field in the vicinity. The 
visualization allows to formulate the following lemma (ana- 
lytically proven in Appendix f). 

Lemma 3 (nearly straight solution, analog of PII-L5,6): 
For the straight line string the singular points of the vector 
field q 3 ±(e3) are: at s = 1 two nodes on the equator (fig'E 1 ), 
at s = 2 saddles on the northern and southern poles and 
four nodes on the equator (fig^ 3 ), at s > 2 multi-saddles 
on the northern and southern poles and 2s nodes on the 
equator (s = 3, fig^O- During the evolution the nodes 
move along the equator, while the singular points in the poles 
stay fixed and saddle patterns rotate around them. After a 
small deformation of the string from straight configuration 
the nodes move along a common trajectory in the vicinity of 
the equator; at s = 2 the saddles move in small loops near 
the poles; at s > 2 the multi-saddles are unfolded to (s — f ) 
non-degenerate saddles moving near the poles. After the 
lapse of time Ar s = n/s = (period of evolution) /2s the vector 
field reverses its direction and the pattern of singularities 
returns to the initial state. During this time the equatorial 
singularities move to the neighbor ones, pole singularities at 
s = 2 perform one revolution along the loops. 

The discussion of the properties of Gribov's copies in 
quantum theory can be found in Appendix 3. 

2 Arbitrary spin case 

2.1 Classical mechanics 

Let's define 7 = {L {2) y 1/2 . 

Lemma 4 (structural, analog of PII-L2): 02 expansion in 
the vicinity of straight-line string has a form 

d 2 

1 — n~> ( 2 ) 

|ai| 4 " v ' 

n> 1 

where P n are polynomials of oi, E_ , S-, d, f, g' , their con- 
jugates and 7,7 _1 . The polynomials possess the following 
properties: (i) each monomial in P n contains, counting pow- 
ers, n variables from the group (E_, S—,d, f,g' and their 
conjugates) and 2(n — 1) variables from the group (oi,a*); 
(ii) each monomial contains at least one variable from the 
group (E±,S±); (iii) variables a± and enter in even pow- 
ers only; (iv) each monomial, containing S-, also contains 
{at) 2k with k > 1; (v) P n = 0(7~ [(n+1)/21 ) at 7 -> 0. 

For the purposes of further consideration it is convenient 
to extract from 02 a common phase multiplier a 2 /|ai| 2 and 
to define 02 = a.2 ■ a 2 /|ai| 2 , so that \a%\ = 1 0:2 1 and 

Q2 = ElJ^' C^+W. (3) 

n>l 

Explicit expressions for the first three polynomials P„ are: 



Pi = -E+ + S+/ 7 , 

P 2 = af f S_ + ai 2 g>* £+ - af f* S_/ 7 - g" S+/y, 

P3 = -af f* g' £_ - ai 2 a* 2 /* g'* S_ - af d* E_ 2 
-ai 2 a* 2 //* S + - ai 4 9 '* 2 S+ - ai 2 a* 2 d£_ S+/2 
+af /* g> 5-/7 + <H 2 af f* g" S_/ 7 + 2a* 4 d* S_ S_/ 7 
+ai 2 af d£ + S_/2 7 - af d* S_ 2 / 7 2 + ai 2 af f f* S+/ 7 
+ai 4 g" 2 S+h + ai 2 af cZE_ S+/2 7 + ai 2 af 7 S_ S+ S+ 
-ai 2 afdS- S+/2f - ai 2 af S + S- S+ - ai 2 af S_ S+ 2 
+ ai 2 af 5_S+ 2 / 7 . 

2.2 Quantum mechanics 

Spin operator was defined in Part II either by components 
in external reference frame (upper index, S 1 ) or by compo- 
nents in the special coordinate system e*;, related with string 
dynamics (lower index, Si). The commutation relations: 

[S\S j ] = ieij k S k , [Si,Sj] = —ieij k Sk, [S\Sj] = 
correspond to raising/lowering operators 

S ± = S 1 ±iS 2 , [S 3 ,S ± ]=±S ± , [S + ,S~} = 2S 3 , 
S± = S 1 ±iS 2 , [S a ,S±] = =F5±, [S+,S-] = -2S 3 , 

i.e. S raises S 3 by 1, S~ lowers S 3 by 1; S+ lowers S3 
by 1, 5*- raises S3 by 1. In further mechanics only low- 
index operators Si will participate. In Part II the oscillator 
variables were composed to the operators E± , ^3 with com- 
mutation relations [A^ 2 \ E±] = ±E±, i.e. E+ raises A 3 by 
(2) 

1, E_ lowers A 3 by 1. These operators define a constraint 

(2) 

X3 = S3 — A3 , representing a symmetry of the mechanics, 
under which the state vectors are invariant: X3|*) = 0. We 
see that the operators E± and S± have similar commutation 
relations with X3- [X3,E±] = TE±, [xs,^] = TS±, so that 
a linear combination of E+ and S+ lowers X3 by 1, while a 
linear combination of E_ and S- raises \3 by 1. We will 
also say that E+ and S+ have X3-charge (—1), while E_ and 
S- have X3-charge 1. Similarly the concept of charges, in- 
troduced in Part II, can be extended to other operators from 
oscillator and spin parts of the mechanics. 

Matrix elements of spin components do not depend on the 
representation of the algebra and have well known form |3]: 



(S(S 3 - 1)5' 3 |5' + |S5'35' 3 ) = y/s(S + l)-S 3 (S 3 -l), 
(S(S3 + 1)S 3 \S-\SS3S 3 ) = y/S(S + l)-S 3 (S 3 + l), 

all other elements vanish. Concrete representation of quan- 
tum top is described by Wigner's functions 53 IS ED- 
Here S characterizes the eigenvalue of Casimir operator 
S'S 1 — SiSi = S(S + 1), commuting with all spin compo- 
nents. S is integer for single-valued representation of SO(3) 
and half-integer for double-valued one. In Appendix 2 we 
show that in our problem only integer values of S should be 
taken. 

Ordering rules. We use the same ordering of elementary 
operators as defined in Part II, Table 1. In the products 
of spin components S± , commuting with elementary opera- 
tors, we select the ordering S+S-, so that S3-raising oper- 
ator S- stands on the right and annulates the states with 
maximal spin projection S3 = S. We remind that classi- 
cally 53 = S — P 2 /2-7T corresponds to one- modal solution, 



straight line string, associated with the leading Regge tra- 
jectory. Quantum analog of © is constructed as follows: 



A' 



02 = n 



~-2n + l . p 
1 • J n 



(4) 



where hi = afai+ci. Here we introduce a constant term ci, 
whose contribution vanishes on classical level (in the limit of 
large occupation numbers afai), and add analogous terms in 
quantum definition 7 = (L^ +C2) -1 / 2 and in the definition 
of mass-squared operator, which we fix as follows: 

p2 



= L {2) + 2a 2 a+ 



■ c 3 . 



(5) 



Note that this definition differs from one used in Part II 
for S = case. Later we will show that this definition is 
preferable for the description of arbitrary spin case. 

Lemma 5 (convergence, analog of PII-L7): 



{h {2) =Ni,S\ :P„ 



'(2) 



N a ,S) = 0, if 



n > min{l + (Ni + N 2 )/2, (4 + 2(JVi + N 2 ) + 45)/5}. 

The resulting spectrum (P 2 /2iv y S) is shown on figO] 
The spectrum has common features with the upper part 
of (Lq 2 ' , A^) spectrum, shown on fig|21 The beginning 
of the spectrum (P /2tv,S) consists of three almost linear 
Regge trajectories. There is a 2-unit gap between the first 
and the second trajectories. The third trajectory starts 
at 5* = 1 level. For the next trajectories the degener- 
ate states of figEl become splitted on figEl The states at 
(P 2 /2n, S) = (3, 1) and (4,0) comprise two numerically close 
pairs with P 2 /2tt = 3,3.0046 and P 2 /2tt = 4,4.0066, while 
the other states on figEl are non-degenerate (not counting 
trivial degeneration for the upper- index S 3 = -S...S and 
direction of P^). The spectrum is computed for the values 
ci = 2, C2 = 4, C3 = 0. Smaller values of ci,C2 correspond 
to higher non-linearities in the spectrum, while larger values 
of ci , C2 make the spectrum more linear and closer to the 
spectrum of (L^,A^ > 0). This behavior is characterized 
by the following lemmas, where for clarity of statements we 
fix C3 = 0. 

Table 1: eigenvectors with P 2 /2-k £ Z. 



P 2 /2tt = L (2) 


S = S3 = A (2) 
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Lemma 6 (analog of PII-L9): the states from Tableware 
annulated by a 2 and have integer- valued P 2 /2n. 

Lemma 7 (leading term): let's keep in only the leading 
(l/ni)-term: |n=l = (~ + S-f^/hi- The states from 
the first two Regge-trajectories are annulated by a? |n=i- 
In this approximation the first two Regge-trajectories have 
integer-valued P 2 /2n and are linear: P 2 /2n = S + k, k — 
0,2. 



Lemma 8 (heavy vacuum modes): in the limit 

1 « c\ « c 2 « c\ the spectrum of (P 2 /2ty, S) tends 
to the spectrum of (L^,^ 2 ' > 0). 
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Fig.2. Spectrum (L (2) ,A (2) ) 



Fig. 3. Spectrum {P 2 /2n, S). 
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Remark (other ordering rules). We have fixed the ordering 
a.20t2 in the definition of mass-square operator after test- 
ing a number of alternative rules: 0^0,2, 0^02, etc. The 
alternative rules produce less regular spectrum (P 2 /2n, S). 
Particularly, for the ordering 0^02 the operator S+, present 
in a 2 ~ (— £+ + S+ /j), does not annulate the states on the 
leading Regge trajectory S3 — S = P 2 /2tv, and introduces 
additional non-linearity to this trajectory. For the ordering 

0,20,2 the operator 01, standing on the right in 02, annulates 

1 (2) \ 

the vacuum state at each S-level: a% \Lq — 0, S) =0, there- 
fore P 2 \L^ — 0, =0 for any S, the resulting spectrum 
looses Regge behaviour. In other tested definition we per- 
formed a search of the pairs a*ai in P n on the classical level, 
and later replace them a^ai — > hi. This definition creates 
almost linear Regge trajectories for small S, possessing more 
non-linearity at large S. 



Conclusion 

In this series of papers we have demonstrated a possibil- 
ity of construction of the quantum theory of Nambu-Goto 
string in the space-time of dimension d — 4. The general 
approach is the selection of the light-cone gauge with the 
gauge axis related in Lorentz-invariant way with the world 
sheet. In this approach the Lorentz group transforms the 
world sheet together with the gauge axis and is not followed 
by reparametrizations. As a result, the theory becomes free 
of anomalies in Lorentz group and in the group of internal 
symmetries of the system. The constructed quantum theory 
possesses spin-mass spectrum with Regge-like behaviour. 

Certain problems are still present in this theory, which 
however do not hinder its implementations, e.g. for the con- 
struction of string models of the hadrons. The results, pro- 
duced by the theory, are influenced by ordering of operators 
and other details of quantization procedure. The theory does 
not contain algebraic anomalies, but possesses the features, 
which can be called spectral anomalies. Particularly, Hamil- 
tonian P 2 /2w, classically generating 27r-periodic evolution, 
in quantum theory is influenced by ordering rules and does 
not have strictly equidistant spectrum. This fact does not 
create problems for the hadronic models, where this spec- 
trum is subjected to phenomenological corrections and ex- 
perimentally is not strictly equidistant as well. The theory 
also possesses a topological defect, appearing as a discrete 
gauge symmetry, identifying the points in the phase space 
(Gribov's copies). This classical symmetry is related with 
discrete non-linear reparametrizations of the world sheet and 
is not preserved on the quantum level. In our construction 
we use the expansion series in the vicinity of one Gribov's 
copy, by these means distinguishing it in the quantum the- 
ory. We have also shown that the leading term of the expan- 
sion, which has a minimal ordering ambiguity and is easier 
for computation, is sufficient to reproduce Regge behaviour 
of the spectrum. Therefore, practically one can keep this 
term to describe the main effect and include further terms 
in the form of phenomenological corrections, together with 
the contributions of other nature [16-21]: gluonic tube thick- 
ness, quark masses, charges, spin-orbital interaction, etc. 
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Appendix 1: proofs for the lemmas. 

Q In the space of variables M = (e^, q n ±, Pn±), introduced 
in Part II, the subgroup generated by X3 ac ts trivially, and 
the whole gauge group of Xi-constraints becomes equivalent 
to SO(3)/SO(2) ~ S 2 . This group defines the vector fields 
qn±{ez),Pn±{ei.), tangent to the sphere |e3 1 = 1. Each of 
these vector fields defines an embedding / : S 2 — > TM(S 2 ) 
of the sphere to its tangent bundle. Let G be an image of 
this embedding for the vector field <f s j_. It is a smooth 2- 
dimensional submanifold of 4-dimensional tangent bundle. 
The gauge condition q s ±_ — defines another 2-dimensional 
submanifold F, which can be identified with original sphere 
S 2 . Due to a theorem p. 525, the index of intersection for 
these manifolds coincides with the index of singular points 
for the corresponding vector field. Then we trivially lift this 
construction to the whole space M. 

H The orbits of Gribov's copies are defined as zero level 
curves of the function F(e- A ) = (q n x p n , es) ~ i{a n x a* , ez). 
This function does not depend on phase rotations, thus the 
orbits are not changed after the multiplication by the phase 
factor. It's easy to verify that in both definitions the evolu- 
tion a n — > a n e~ lnT is equivalent to the motion of the marked 
point O along the supporting curve Q(o) under the action of 
uniform shifts: a — * a+r. The only difference is that for PII- 
Eq.(12) these shifts are performed in lcg-parametrization, 
while (Q corresponds to the natural parametrization. 

[3] The structure of singularities near the poles has been 
described in PII-L5,6. On the whole sphere of ea for 
the straight string solution Q{a) = (cos a, sin a, 0) the in- 
tegral Q can be evaluated in terms of Bessel functions: 

a n = H)"yV8(J n+1 (np)e l(n+1)Q + J^np)^^- 1 ^, 
-iJ„+i(np)e l(n+1)a + iJ„_i(np)e !(n - 1)Q ,0), where e 3 = 



(pcosa,psina,p z ), < p < 1, p 2 + p 2 = 
1, n = 1,2... The real part q n = Re(a n ) = 

yj 7r/8(J„+i(np) cos + + J„_i(np) cos <j>-, J n+ i(np) sin — 
J n _i(np) sin0_, 0), where 0± = (n ± l)a — nn/2. Ze- 
ros of g„x = q„ — (q n ez)ez correspond to q n x e 3 = 
(g^ejj, —q^e^q^e 2 — <& e \) = 0. The solutions disjoin to the 
following branches: 

(bl) e\ = 0, equator: q\e\ - q 2 n e l A = 0; 

(b2) e| ^Q,qi = ql = 0. 

For the case (bl) p = 1 and after simplifications we have 
( J n -i i n ) — Jn+i (n))-sm(na— nn/2) — 0, i.e. a = n/2+nk/n, 
k — 0...2n — 1. Considering the expansion of the vector field 
e.g. near the point a = n/2 in the local coordinates (Aa, p z ), 
we have g nX ~ (-n(J„_i(n) - J n+ i(n))Aa, -(J n -i(n) + 
J„+i(n))p z ). Because of the property J n -i(n) > J„+i(n) > 
(see below) we conclude that this singular point is the 
node. For the vector fields, corresponding to the straight 
line string, the evolution is equivalent to the rotation about 
z-axis. The evolution during Ar n = n/n matches the node 
on the equator to the neighbour one. Let's consider small 
perturbations of the string from the straight configuration. 
Non-degenerate points persist and preserve their type after 
the deformation. From P1I-L3 we know that singular points 
move along zero-level curves of the function F(es) = (q n x 
p n , e%). For the straight line string this function is 



F ■ 



p cos a 
„i 



p sin a p z 

a\ 
,2 



"7, 







pz{J n ~i{np) - J n+1 (np)). 



Near the equator F ~ p z , i.e. the equator is non-degenerate 
zero level curve of this function: dF/dp z 7^ 0. Parameter- 
izing zero level curve by the longitude: 63(a), after small 
deformations: F(ez + Se-j) + SF(ez) = 0, for the shifts of 
the latitude we have: Sp z = -(dF/dp^SF. Therefore, 
the curve will be slightly deformed in the vicinity of the 
equator. It was proven in the Part II, that the vector field 
reverses its sign after the lapse of time Ar n . During this 
time the pattern of singularities returns to the initial state. 
Thus, the nodes near the equator should move to the neigh- 
bor ones, to satisfy the continuous limit with the straight 
string configuration. 

For the case (b2) p < 1 and we have a system 



J n -i(np) aax(<t>- 

Jn + l{np) + J n - 



- + <t>-)=Q, 
i(np) cos(4> + + 1 



0. 



Here we again have two branches: 

(b21) J n+1 (np) = J n -i(np) = 0; 

(b22) 0+ = Trfc, J n+ i{np) + (-) fc J n _i(np) = 0. 

In the case (b21) we have only solution p — at n > 1. 
It gives pole singularities described in PII-L5,6. Near this 
solution J n +i « J n -i « 1 and in the passage around this 
point the vector (g^,^) ~ p n-1 (cos0_, — sin0_) performs 
(n — 1) revolutions in the direction opposite to the pas- 
sage. It corresponds to the saddle for n — 2 and multi- 
saddle for n > 2. In the complex plane of z = pe la the 
mapping z — > q — q„ + iq\ near this point can be rep- 
resented as q ~ (z*) n_1 . Considering small deformations of 
the string from straight configuration, we introduce a correc- 
tion (z*)™" 1 + 8q, where Sq(z, z*) = 5qo + Sqo z z + Sqo z * z* + ... 
is small together with its derivatives, and in small vicinity 
of the singular point we keep only the constant term Sqo- 
At a given point there is 1-1 correspondence between the set 
(q n ,Pn)o and the shape of supporting curve, and the arbi- 
trary small variations 8qo can be reproduced by variations 



of supporting curve. That's why we can apply the same per- 
turbation analysis as for general vector fields, in spite of the 
fact that we consider the fields of special form. New singu- 
lar points are z k = (_«5^) 1 /(«-i) e 2^/(n-i) > fc = g...n _ 2 . 
Expansion of the field in their vicinity has a form q = (zj* + 
Az*) n ~ 1 +Sq ~ (n-l)(4) n_2 Az*, therefore we have (n-1) 
saddles. Considering their evolution: q — * q cos nr+psin nr, 
where p ~ i(2*) n_1 + 5p, we have (z*) n ~V nT + Sqcos nr + 
Spshmr = 0, i.e. z k = {-8z\ - fc 2 *e 2 " lT ) 1/(n - V e *«ik/(n-i) ) 
where 5z± = (5q — i8p)/2, 8z2 = (5q + iSp)/2. If n = 2, 
we have a single saddle moving along a closed loop with a 
period Ar 2 = n/2, as described in PII-L6. If n > 2, the 
dynamics is more complex, see fig'0] Particularly, if 8z\ = 0, 
8z2 7^ 0, 71 — I saddles move along a common circle, match- 
ing the neighbour one after the lapse of time Ar n = n/n. 
If Szi 7^ 0, Sz2 — 0, the saddles stay fixed. Generally, if 
< \5z2\ < \Szi\ « 1, the origin z = is placed outside the 
circle Sz\ + Sz2e~ 2lnT , and the complex root of (n — l)-th 
order gives (n— 1) saddles moving along disjoint loops with a 
period Ax„. If < \8z\\ < \St2\ « 1, the origin is inside the 
circle and the saddles move along a common loop. The sep- 
aratriss \8zi\ = |<5z2j for n > 3 is structurally unstable and 
under the influence of higher order corrections is unfolded 
to further complex cases. 




Fig. 4. Unfolding of multi-saddle, 



6. 



Note: this phenomenon resembles the bifurcation of de- 
generate critical point of the function / = x 3 — 3xy 2 |S] 
(known also as "monkey saddle" |7|), and the structure of 
phase flows near high order resonanses of Hamiltonian sys- 
tems mi- 

In the case (b22) we have a = nk/2n, k — 0...4n — 1, 
J n+ i(np) — J n - 1 (np) — for odd k and J n+1 (np) + 
J n -\(np) = for even k. The solution p = 0, n > 1 has been 
already considered. There are no solutions in the interval 
< p < 1 due to the inequality J n _i(np) > J n +i(np) > 0, 
which follows from easily provable relation for Bessel func- 
tions: J n -i{x) > J n (x) > J n +i(x) > ... > for < x < n, 
n= 1,2,3... 

|1] The equations (6) of Part II in the arbitrary spin case 
can be rewritten as 
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13 = a 2 a 2l 2 - (3 2 /2, 

where ri\ = a^ai, z = E+ + /*a2 + g'*a2 + d'a^ 2 + da2a 2 /2, 
■ = \J P 2 /2tt. Substituting the expansions 



1 



2 



E 

n>l 
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/3=> zSr 
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we obtain the following recurrent relations: 



Pi = -S+ + — , 6i=0, 



Km<n-2 



±dal 2 al ^ P n _ m _iP^ + —b n -i, n > 1 



: a* 2 ai7 2 ^ P*_ m Pm - 5 



Km<?i — 1 



2 / ,, 

Km<n — 1 



in-ml>m, n > 1. 



Then the properties (i)-(v) can be easily proven by induc- 
tion. The property (i) also follows from the scaling symme- 
try of the system (JSJ: formal transformations {ni — > \f(J-n\, 
v — > C ■ v, v G (d, /,<?', E_, S 1 - and their conjugates)} and 
{ni — > C • ni, ai — > C • cii} preserve 0,2, (i\ the property (in) 
follows from reflection symmetry: a\ — > — 01 preserves 0,2, /3. 
The property (ii) follows from the fact that a2 = /3 = is 
the solution of (|SJ at E± = S±. 

Let's consider a monomial of P n , substitute the def- 
initions of elementary operators d = d+ — d_, / = 
/+ ~ /-, 9' = 9- - 9+ (see PII-Eq.(7)), expand it 
to a number of (secondary) monomials and normally or- 
der each one. Let n(df),n(ff),n{gi),n(aif) be the num- 
(2) 

bers of L -raising operators in the secondary monomial, 
n(di),n(fi),n(gi),n(aii) be the numbers of Lg -lowering 
operators. Using the definitions of the charges (Table 1 
in Part II), we see that the secondary monomial is repre- 
sented as a diagonal with offsets A^L^ = 4n(dt)+2n(/f) + 
2n(g' t )+n(a n ), A,L, (2) 



H ±. = 4n(d l ) + 2n(f l ) + 2n(g{) + n(a n ) 



<2) 



in the matrix, displayed at fig^: = ATi|mon|Lg 

N 2 ) ~ (l[ 2) = Ai-A T L 2) |dia 5 |L 2) = Ni-AiL^). Here 

(2~) 

diag represents L -neutral part of the monomial, consisting 
of the operators £± , S± , 7, and non-zero value of this matrix 



element corresponds to Ni ~ A^L 



(2) 



iV 2 



A|4 2) > 0. 



N m N m +2Aj; 




(a) (b) (c) 

Fig. 5. To the proof of lemma l5l 



Then, using the inequality \A$ \ < and charge prop- 
erties of E±, see figEJs, we conclude that non-zero value for 
the considered matrix element corresponds to the following 
restrictions on the number of E±-operators: n(E+) < 2(N2 — 



"(2) 



Ai4 2) ), n(E_) < 2(N 2 - AiL£>), i.e. N 2 > + A E , 

Ae = max(n(E+), n(E_))/2. As a result, non-zero entries 
in the diagonal, shown on figE^, receive an additional offset. 
After summation over all secondary monomials non-zero en- 
tries of the matrix are localized in the region, shown on figHt. 
HereiV m = A T L 2) + A^ 2 ' = 4n(d) + 2n(f) + 2n{g)+n( ai ) 
and n{x) represents the number of variables (x,x*) in the 
primary monomial. Particularly, non-zero entries are local- 
ized at Ni + N 2 > N m + 2A E . For the number of Si- 
operators we have the following restrictions: n(S+) < 2S, 



(2K 



-(2) 



n(S-) < 2S. Then, using the property (i) from the struc- 
tural lemma: n(S±, £±, d, f, g) = n, n(ai) = 2(n — 1), we 
construct the following system of 11 linear inequalities and 
1 linear equation: 

4n(d) + 2n(f) + 2n(g) + n(S + ) < Ni + N 2 + 2 - 2n, 
4n(d) + 2n(/) + 2n{g) + n(S_) < iVi + N 2 + 2 - 2n, 
0<n(S+)<2S, 0<n{S-)<2S, 
n(S+) > 0, n(E-) > 0, n(d) > 0, n[f) > 0, n(g) > 0, 
n(d) + n(f) + n(g) + n(E+) + n(£_) + n(S+) + n(S-) = n. 

The inequalities define a convex polyhedron in 9- 
dimensional space (n,Ni + N2 + 2, S, n(d), n(f) + 
n(g),n(Y,+ ),n(Y,-),n(S+),n(S-)). The linear equation de- 
fines its 8-dimensional slice. Projecting it to a 3-dimensional 
subspace of parameters (n,Ni + N2 + 2, S), we obtain a 
convex polyhedron, which represents the region of existence 
of the solutions. Solving the system by Mathematica, we 
have n < min{l + (iVi + N 2 )/2, (4 + 2(N 1 + N 2 ) + 4S)/5}. 
This condition defines a region, where non-zero entries of 
the matrix element considered in the lemma are located. 

El The states, annulated by ajj", have P 2 /2n — G Z, 
at C3 = 0. Due to the property (ii) of structural lemma 
each monomial in P„ contains at least one operator from 
the group (E±,S±). After the normal ordering some num- 

(2) 

ber of L -lowering operators can also appear on the right 
of E±. Due to the charge properties and the property (iii) 
of lemma the Lq -lowering operators together decrease 



'(2) 



(2) 



2k, k G 

(2) 4 (2K 



Z, k > 1. The vacuum state 
(0, 0) is annulated by all these 
(1,1) is 



(2) 4 (2K 



(P 2 /2tv,S) = (L^>,A 

operators. The state (P 2 /2ir,S) = (L^ 1 , A 
annulated by E±-operators, together with other states on 

the leading Regge trajectory, see figEl If there are any L^- 

(2) 

lowering operators at the right, they decrease L by > 2 
and annulate this state. It is also annulated by S- and 
generally is not annulated by S+. However, due to the prop- 
erty (iv) in lemma 0] each S+ enters to a\ together with 
a 2 , which annulates this state. For the state (P 2 /2iv, S) = 
(Lq 2 , A^) = (2,0) the terms with S± are inactive. E± an- 
nulate this state, see fig|21 If there are any Lq -lowering 
operators at the right, they can transform this state only 
to vacuum, and the vacuum is also annulated by E±. Sim- 
ilarly, on the states (P 2 /2tt,S) = (Z^, 2) ,^ 2) ) = (4,0) the 
terms with S± are inactive. There are 3 states on the level 



(Z/q , j4g ) = (4,0), 1 state on the level (L^>,A^>) = (4, 1) 
and 1 state on the level (Z^ 2) , Af ] ) = (4,-1), see fig|U As 
a result, the operator E+ annulates 2-dimensional subspace 
of (L^\ A^) = (4, 0) and one state in this subspace is also 
annulated by E_. Using the definition of E±, it's easy to 
verify that E±|2i2_i) = 0. If there are any -lowering 

operators at the right of E± , they can transform this state to 

( 2) 

the states with L = 0, 2, which are also annulated by E±, 
see fig|U There are 2 states on the level {L^\ A^) — (3, 1): 
I3), J 2i 1_ 1 ^ . Both states are annulated by E±- If there 
are any Lq 2 ' -lowering operators at the right, they transform 



these states to L Q = 1, also annulated by E±. Both states 
are annulated by S-. Generally they are not annulated by 
S+, and due to the property (iv) in lemma^Jeach S+ enters 
to q.2 together with a\, which annulates one of these states: 

|l3>. 



In a^j n=1 = (—Yl—+S-/ r f)/hi the operator £_ decreases 
A$ by 1, the operator S- increases 5*3 by 1. Let's consider 
the states A 2 = S3 = S > 0, annulated by 5*-. On fig|21we 
see that £_ annulates the leading Regge trajectory, because 
there are no states immediately below this trajectory. Con- 
cerning the second trajectory < A^ — — 2, one state 
at (L f (2) , A (2 } ) = (2, 0) and two states at {L {2) , Af } ) = (3, 1) 
are annulated by £_ due to the same reason, while the se- 
quence of double states 2 < A^ = — 2 is transformed 
by the operator £_ to the sequence of non-degenerate states 

1 < A^ — Z/q 2 ' — 3, thus annulating one state from each 
pair. As a result, the described states are annulated by 
a+| n=1 and in the considering approximation these states 

have P 2 /2ir = from here we obtain the statement of 

the lemma. 

1 4 = £P+/n 2n -\ where P n = 0( 7 - [(n+1)/21 ) 
due to the property (v) of lemma 0] Using the 

definitions h\ = a^a\ + ci, 7 = (L^ + C2)" 1 '' 2 , 
in the limit 1 « c 2 « C2 « cf we have ot%\ n ~ 
4 n+1)/2]/2 /c?-\ i.e. a+| n=1 ~ c\ /2 /c,»l, a+|„ =2 ~ 
C2 /2 / C f«l/ Cl «l, a+\ n > 3 ~ cr+^/c 2 - 1 < 
4™ +1)/4 /c 2n - 1 « l/c'j 1 " 2 « 1. Writing explicitly the lead- 
ing contribution: a\ — S- ■ c^ 2 /c\ + o(l), we see that the 

states with S = S3 = A^ are annulated by in the de- 
scribed approximation, while other states receive large con- 
tribution to P 2 /2n and are shifted away from the origin 
of the spectrum. Due to this property, the spectrum of 
(P 2 /2tt, S) tends to the spectrum of (L ( g ) ,A^ ) > 0), pro- 
vided that C3 = 0. 

Appendix 2: super- select ion rule. 

In several earlier works the quantum top appeared in the 
context of relativistic string dynamics. Particularly, the pa- 
per Hj] describes a special family of motions in a theory of 
Y-shaped strings, globally isomorphic to the phase space 
of the top R 3 x SO(3). In the paper 10 a similar fam- 
ily was selected in the phase space of open strings. In the 
present paper we also have SO(3)-factor in the topology of 
the phase space. As a manifold, it has a fundamental group 
iri(SO(3)) = Z2, i.e. it possesses two homotopically non- 
equivalent classes of contours. The contour, representing 
one complete revolution, cannot be contracted to a point 
(unit of the group), while two revolutions are contractible 
|14) . This fact gives a possibility to construct ray repre- 
sentations of rotation group, corresponding to integer and 
half-integer values of spin [15) . The existence of half-integer 
spin and associated double-valued representations of 5*0(3) 
is also explained in terms of path integration technique: non- 
homotopic contours can contribute to path integrals with 
different phase factors [131 15]. 

The question on the availability of half-integer spin is 
subtle and requires special consideration in each particu- 
lar problem. Now we will consider a model example, rep- 
resenting the main features of our problem: the mechanics 
of 3 particles in 3-dimensional space, (x n ,p„) G R 9 x R 9 . 
Let's apply the following canonical transformations: Q = 
dp„ A dx„ = dP A dX + |d(5 x e*) A det + dp ni A dq ni , 
where X = ^m n £ n /^m n , m n are masses of the parti- 
cles, P = V^pn, q n = x n — X, p n = p n — m n /m3 ■ ps, 
S = Y]qn x p n — ^2<ln x Pn, Si is an orthonormal basis, 
pni = pn&i, q n i = QnSi, and everywhere the sum over re- 
peated indices is assumed. Let's fix the basis with respect 



to the particles. The following cases should be considered: 

(si) <fi and 52 are linearly independent. In this case we can 
select e.g. ei TT <?i, e 2 TT <72± = 92 — (g*2gi)/<7i 2 • qi- The last 
term in the symplectic form can be written as dpu A dgn + 
dp2iAdq2i+dp22/\dq22. Herein > 0, 922 > 0. The obtained 
phase space has topology (R 9 ) x (R 6 x 50(3)), where the first 
factor represents spin-momentum space, the second factor is 
configuration space, 5*0(3) is the configuration space of the 
top. 

(•S2) q\ / 0, q 2 / 0, q\ TT <7*2- Let's direct e\ TT <7i- The 
symplectic form can be written as fl — dP A dX + d(S x 
ei) A dei + dpu A dqu + dp2\ A dq2\. Here Set = 0, gn > 0, 
Q2i > 0. The topology of the phase space is R 10 x TM(S 2 ), 
the configuration space is R 5 x S 2 . Here the tangent bundle 
TM(S 2 ) represents the phase space of rotator. 

(53) 9l 7*= 0, q2 7*= 0, qi \{ 9*2- Corresponds to the same 
mechanics as (82), but 921 < 0. 

(54) q\ 7*= 0, 52 = 0. The same mechanics as (sa), but 
Q2i = and P21 omitted: phase space R 8 x TM(S 2 ), config- 
uration space R 4 x S 2 . 

(ss) 91 = 0, q\ 7^ 0. The same mechanics as (S4), but 
g22 = and P22 omitted. 

(fi6) 9i = 92 = 0. In this case S = and 57 = dP A dX. 
The phase space is R 3 x R 3 , scalar particle. 

The archetypes (top, rotator, scalar) are essentially the 
same as described in Sec. 2 of Ref.|13|. The important dif- 
ference is that in our example they appear not as isolated 
cases, but as a stratification of the configuration space: 
R 6 (9i,g*2) = uf =1 Si. The strata S2-6 have zero measure and 
codim{s2,'i) = 2, codim{sA,5) = 3, codim(se,) — 6 in R 6 . The 
exclusion of these zero measure subsets changes the topolog- 
ical structure of configuration space. Note that the original 
configuration space is simply connected and quantization of 
the problem in original representation gives only integer val- 
ues of spin. All closed loops in the original configuration 
space are contractible. The contours in this space can be re- 
lated by a homotopy, intersecting the strata S2-6, as a result, 
all the contours should contribute to path integrals with the 
same phase factors. Therefore, only single- valued represen- 
tation of 50(3) can be used. We conclude that half-integer 
spin in the considered model example appears as an artefact 
of the representation and should be rejected. 

In our problem constraint S3 = A3 £ Z selects only inte- 
ger 5, however, with formal modifications A3 — > A3 + 1/2, 
classically vanishing after the recovery of Planck's constant, 
it's possible to introduce half-integer spin to the theory. The 
above argumentation favors against it. During our construc- 
tion we reject zero measure subsets of singular cases (see 
Part II). The original configuration space of open string the- 
ory is simply connected, all contours in it are homotopic to 
each other. Therefore we should not introduce half-integer 
spin in this theory. 

In the problems EJ3 submanifold R 3 x 5*0(3), rep- 
resenting particular type of string motion, was selected in 
the original phase space, representing all possible motions. 
The original configuration space possessed trivial topology: 
x a {cr) G O°°[0, 7r] for open string and x^(cr) G O°°[0, n] with 
x li( n ) = = f° r Y-shaped string. Therefore, half- 

integer spin representations have right to exist in frames of 
the restricted models [jU I1U| . but would disappear in full 
quantum theories. 

Appendix 3: Gribov's copies in quantum mechanics. 

The Gribov's copies of string theory can be separated to 
those which belong to the same xo-orbit and are transformed 
one to the other by the quarter period of the evolution, and 



those which belong to different orbits. The factorization of 
the phase space with respect to xo-transformation identifies 
the copies on the same xo-orbit. 

The discrete symmetry D4 — exp[— i(7t/2)Lq 2 ], mapping 

the Gribov's copies on the xo-orbit one to the other, is as- 

(2) 

sociated with the quantum number Q — L mod 4. In 
the whole phase space, containing the external variables 
(X, P), xo-orbit is not closed, but is periodically repeated 
as X — > X + 2Pn, n G Z. Considering the evolution in 
this space, the generator of discrete symmetry should be ex- 
tended to Di = exp[i(7r/2)(P 2 /27r - L {2) )]. Classically, for 
the Gribov's copies near the poles this evolution coincides 
with D'l = exp[i(7r/2)xo]- Indeed, the evolutions D'4 and 
D'l connect the same Gribov's copies, but along different 
trajectories, particularly, D'4 does not contain Si -variables, 
influencing the direction of gauge axis, and D'l does. For the 
Gribov's copies near the equator the evolutions D'4 and D'l 

connect not the same, but neighbor copies. These formal dif- 

(2) 

ferences become crucial in quantum mechanics, where L 
has integer- valued spectrum and M 2 /2% = -\-2ctia\ has 
not. Note that the number Q is associated with L and D4- 
symmetry, defined in the phase space of internal variables, 
independent of (X, P). 

Let's write xo = (P 2 — M 2 )/2tv, where M 2 is a function 
of the internal variables and has a discrete spectrum. Let's 
consider the Hamiltonians of interaction H lnt , which gener- 
ate transitions in this spectrum. Generally H{ n t should not 
commute with xo: this would mean that H[ nt is constant on 
Xo-orbit. Also it should not commute with D'4, otherwise it 
would have the periodic properties, which select too narrow 
class of Hamiltonians. It was shown in paper |20) that not 
Hamiltonian itself enters to the physical vertex operators, 
but its average over a period of evolution. This averaging 
actually inserts a projector onto the subspace xo = 0. In 
a similar way we can construct an operator J, commuting 
with xo or D'4. Let's consider two cases: 

(Ja) J does not contain X. In this case [J, xo] = im- 
plies [J, M 2 ] = 0, i.e. J cannot create transitions between 
the states with different masses. Also, [J, D'4] — implies 
[J, Da] = 0, i.e. J cannot create transitions between the 
states with different Q. 

(Jb) J contains X. In this case the operator with [J, xo] = 
generates simultaneous changes of P 2 and M 2 , preserving 
the mass shell condition. In a similar way, J creates transi- 
tions between Q-sectors, even if the condition [J, D'4] — is 
satisfied. 

Note that the operator M 2 itself belongs to the class (Ja), 
as a result, it is diagonal with respect to the quantum num- 
ber Q. The vertex operators generally belong to the class 
(Jb): they contain X-terms, responsible for the change of 
P 2 . E.g. the interaction with a planar wave |20| has a factor 
e lkx , changing P — *• P — k. For such interactions Q-sectors 
are mixed up, therefore no selection rule for Q is necessary. 

After factorization of xo-transformations the remaining 
discrete gauge symmetry is related with the choice between 
different xo-orbits. For nearly straight string they are repre- 
sented by equivalent solutions in the vicinity of the northern 
pole, the equator, and the southern pole. In the obtained 
spin-mass spectrum the first Regge trajectory corresponds 
to the northern pole solution, further trajectories are cre- 
ated by expansion |J5J in the vicinity of the northern pole 
solution. Classically there are gauge equivalent solutions 
near the southern pole and the equator, which possess the 
same (P 2 ,S). In quantum mechanics we do not see these 
additional solutions in the spectrum. Particularly, the first 



Regge trajectory is non-degenerate. We conclude, that the 
equivalence between these solutions is lost in quantum me- 
chanics. The spectra for additional solutions can be shifted 
to the region of large masses, or the quantum expansions 
Q can even diverge on these solutions. Indeed, the classical 
solutions in the vicinity of the northern pole possess large 
n\ and small nj,, k G Z/{0, 1,±2}, providing the conver- 
gence for (l/ni)-expansion J^J. In quantum mechanics the 
convergence of expansion is supported by the normal or- 
dering of operators and the finite number of occupied modes 
in (Lq 2 ^ < iV)-spaces. For the solutions near the southern 
pole ni — > 0, and the usage of (l/m)-expansions is problem- 
atic. One can use (l/n_i)-series to construct a definition of 
mass shell condition, alternative to 0, however these defi- 
nitions will substantially differ on the quantum level and in 
fact will create two distinct theories. For the solutions on 
the equator infinite number of oscillator modes are excited, 
and for these solutions the convergence of the expansions is 
not guaranteed neither on classical nor on quantum level. 
This argumentation explains why the usage of (l/m)-series 
in the vicinity of the northern pole solution preserves only 
one Gribov's copy in the quantum theory. 




Fig. 6. Example: Gribov's copies in the phase space 
of planar rotator. 

Remark. The problem of Gribov's copies in full generality 
was addressed in the paper This paper has consid- 

ered several systems, possessing the Gribov's copies, and dis- 
cussed various approaches to their quantization. In the sim- 
plest example, planar rotator, the phase space M = S 1 x R 
is formed by a periodical coordinate tp mod 2-7T and spin vari- 
able S G R. There is a constraint E — S 1 of a form 
S — Const, whose gauge group G = S 1 identifies all points 
on the circle and therefore eliminates from the theory all de- 
grees of freedom. According to such system does not 
allow gauge fixing condition F of the form / = 0, where / is 
a smooth real-valued function, globally defined on M, which 
has only non-degenerate zeros. Indeed, non-degenerate zero 
of / divides the circle to the parts with / > and / < 0, 
and there is at least one other point with / = on the circle, 
see figJSJi. However, for construction of the reduced phase 
space it is not necessary for / to be globally defined. For 
the computation of Dirac's brackets it's only needed that / 
has been defined locally, in the vicinity of the intersection 
E n F. Generally F can be an arbitrary smooth subman- 
ifold of M, not necessary representable as zero level of a 
globally defined function, see fig|U>. Moreover, for the tech- 
nique |H| I12| only the intersection Ef~) F should be specified 
and the symplectic form of M can be reduced directly to 
E fl F. In the considered example such class of gauge fix- 
ing conditions allows unique intersection between E and F. 
In this low-dimensional example the surface of constraint E 
coincides with the orbit of gauge group G, generally they dif- 
fer. In certain cases the topology of gauge orbit G and gauge 
fixing condition can lead to multiple intersections. Partic- 
ularly, if both are closed submanifolds of Euclidean space 
(with dim(G) = codim(F)), their intersection index equals 
zero, giving an even number of transversal intersections [Sj- 
The same property holds if one is a closed submanifold and 



the other one is a linear subspace. The set EnF can be con- 
nected (otherwise the connected component can be taken), 
and the orbits G of gauge group, to which E is foliated, can 
intersect F several times. In this case the mechanics pos- 
sesses Gribov's copies and associated discrete gauge symme- 
tries. 

Appendix 4: numerical methods. 
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Fig. 7. S^-extension of the algorithm. 

The necessary extensions of the algorithms developed in 
Part II from S = to S > consist in the introduction of 
the additional block structure, related with the quantum 
number S3 £ [— S, S]. Now we need to compute the matrix 
elements (x3 — 0,Q|a2|x3 — 1>Q)i shown on fig[7|at the 
right. Here the large blocks correspond to S3 number, small 
blocks correspond to A3. The required matrix elements are 
localized in the intersection of rows with — S < S3 = A 3 < S 
and columns with — S < S3 = A 3 + 1 < S. Actually, the 
only extension is required in the algorithm ord, procedure 
inject, see Part II. The results of sparse matrix mul- 
tiplication after injection to (A 3 f, Q|q2 | A 3 i, Q^j blocks 
in the left matrix on ngQ should be further passed to 
the right matrix. However, keeping in this matrix only 
necessary elements, we condense it back to the smaller 
left matrix and perform the actual computations in it. 
Each monomial has a definite AA3-charge and therefore 
contributes to one of the diagonals in this matrix, marked 
on figEl (S3 = A 3 \mon ■ 5 ™( s +> S "( s -> | Sa - n(S_) + 
n(S+), A 3 - n(S-)+n{S+) - l) = (A 3 \mon\A 3 - n(5_) + 

n(S+) - l) ■ (S 3 - n(S_) + n(S+)\S[ 
The last matrix element gives 

Tls 3 —n(S_) + -n(S + )<m<S 3 +7i(S + ) V^O^ + 



i(S_) „n(S 



s: w+) |s 3 ). 

coefficient 



- m(m - 
Here 



the 



ris3<m<S3 + n(S + ) \/ S(S + 1) - m(m + 1) . 

case of empty set of indices (e.g. n(S+) = 0) we define 
J] = 1. If S3 + n(S+) > S or S 3 - n(S-) + n(S+) < -S, 
we have Y[ = 0. Substituting S3 = A3 in these for- 
mulae, we have non-zero coefficients only for the values 
A 3 + n(S+) < S, A 3 ~n(S-)+n(S+) > -S, i.e. A 3 changed 
between A 3m in = max{ — S, — S + n(S-) — n(S+)} and 
A 3m ax — S — n(S+). If A 3m in > A 3ma x , the monomial is 
omitted. Then we initialize the sparse matrix multiplicator 
as projector to (A 3 € [A 3mini A 3max ], mod 4 = Q) 
subspace and perform the rest computation as described 
in Part II. As an additional optimization, in the products 
J"J. op^ we have implemented preliminary computation and 



storage of the powers opf for k = 1. 
this procedure by a factor 1.8. 



.kmaxi, accelerating 




Fig.8. Spectrum (P 2 /2tt, S - en) as a function of n. 




Fig. 9. Spectrum (P 2 /2tt, S - eN) as a function of N . 

The determined spin-mass spectrum as a function of the 
correction order n and cutoff parameter is presented 
on figs 18191 The figure |H| shows superimposed spectra 
(P 2 /2ir,S — en) at fixed N = 11, where e is a small con- 
stant and n is changed from 1 to 6 in the direction from top 
to bottom. By the reasons explained in the lemmas above, 
some of the states remain fixed (vertical) with the change 
of n, while the others are rapidly stabilized with n. On the 
figure the spectra {P 2 /2n,S - eN) are shown for n = 3, 
TV = 7. ..20, where in each Q-sector values /Vmod4 = Q were 
selected. When N increases, the most of the states also show 
tendency to stabilization. The spectrum, presented on figOl 
is computed for n = 3 and N — 17, 18, 19, 20 for 4 values of 
Q superimposed in this spectrum. 

The other statistics of the algorithm is presented in the 
following tables: 

Number of secondary monomials in each order n = 1. 



n 


1 2 


3 


4 


5 


6 


num. of monomials 


2 8 


48 


248 


1066 


3844 



Time and memory requirements of the algorithm 





(n = 


3,5 = 6, Q 


= 0): 




N 


8 


12 


16 


20 


subdim 


85x83 


447x438 


1932x1900 


7288x7183 


req. memory 


2.1Mb 


5.6Mb 


83Mb 


1.1Gb 


comp.timc 


0.08 sec 


1.4 sec 


14 min 


15 hours 



In the last table subdim is the dimension of S = 6, Q = 
block in dense matrices, which at large N mainly defines the 
amount of required memory. The computation is performed 
on HP 2GHz Linux PC. 



This figure "grib.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-th/0306010vl 



